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Abstract 

We estimate the effects of non-perturbative physics on the differential distribu- 
tions of infrared- and collinear-safe e + e~ event shape variables, by extending 
the notion of an infrared-regular effective strong coupling, which accounts for 
the non-perturbative corrections to the mean values of several shape variables, 
to their distributions. This leads to 1/Q power corrections over a range of val- 
ues of the shape variables considered, where Q is the centre-of-mass energy. In 
the case of the thrust variable, the leading correction is simply a shift of the 
distribution, by an amount proportional to 1/Q. We show that this gives an 
excellent description of the data throughout a wide range of T and Q. 
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1 Introduction 



The study of power-suppressed corrections to QCD observables has become a lively field 
of experimental and theoretical investigation. On the experimental side, the estimation 
of power-suppressed hadronization or higher-twist corrections is necessary for the accu- 
rate measurement of the strong coupling a s , a crucial parameter of the Standard Model. 
Theoretical work on divergences of the QCD perturbation series, in particular infrared 
renormalons |3J@], and related attempts to define the running coupling beyond perturba- 
tion theory ||, have also led to renewed interest in power-suppressed contributions. 

Power corrections to hadronic event shapes are of particular interest, because they are 
expected to have a characteristic 1 / Q-dependence on the hard process scale Q (the centre- 
of-mass energy in e + e~ annihilation) |4]-|T0||. Thus, in contrast to the 1/Q 2 corrections 
encountered in most quantities, they are still important even at Q ~ Mz, where most of 
the e + e~ data have been obtained. The theory of l/Q-corrections is much less developed 
than that for higher powers; in particular, the operator product expansion does not apply 
in this case, since there are no relevant operators of the corresponding dimension. 

In refs. @,||]|| a treatment of l/Q-corrections to event shapes was proposed, based on 
the notion of an effective strong coupling, a e «, which is approximately universal but differs 



from the perturbative form in the infrared region (see also Such an approach was 

found to be quite successful in describing the powers and approximate magnitudes of power 
corrections to a wide variety of QCD observables, using the low-energy moments of a e s as 
non-perturbative parameters. The comparison with data on event shapes, however, has 
been limited so far to the mean values only ||5|,[12|. 

In the present paper we apply the same approach to the differential distributions of event 
shape variables. The central idea remains that of ref. |J: the emission of soft gluons is 
assumed to be controlled by an effective coupling a eff , different from the perturbative form 
in the infrared region but small enough for terms of higher order in a e Q to be neglected as a 
first approximation. We combine this idea with the treatment of event shape distributions 



developed in refs. ||13| , |T4|| . It was shown there that some shape variables have the property of 
exponentiation, which allows large logarithms to be resummed to all orders in perturbation 
theory. For such variables, the effective coupling assumption means that the leading non- 
perturbative corrections also exponentiate, implying a specific transformation of the whole 
distribution, or at least its logarithmically enhanced part. This has been pointed out in 
ref. [0. 



We concentrate here, as in ref. ||, on the thrust variable, T \TE\. In this case the leading 
non-perturbative effect over a range of thrust values turns out to be simply a shift in the 
distribution, by an amount proportional to 1/Q, modulo logarithmic Q-dependence. The 
shift is just such that we recover the result of ref. || for the mean value. Remarkably, this 
leads to an excellent description of the data over a wide range of T and Q. The only two 
free parameters are the non-perturbative quantity 

ct (/!i) = — / dqa cS {q) (1.1) 



Hi Jo 

which characterises the behaviour of the effective coupling below some infrared match- 
ing scale /jl u and the perturbative coupling a s (M z ). For these quantities we find values 
consistent with those obtained from other data. 
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After a detailed consideration of the thrust distribution, we discuss the relation between 
our analysis and the "dispersive approach" of ref. 0, and then comment briefly on the 
prospects for extending the method to other shape variable distributions. 



2 Thrust distribution 



It was shown in ref. [[14] that the thrust distribution is given to next-to-leading logarithmic 
accuracy by the expression 



1 da 

o "tot dT 



F{T) 



2m Jc 



J q M 



(2.1) 



where the contour C runs parallel to the imaginary axis, to the right of all singularities of the 
integrand, and J^(Q 2 ) represents the Laplace transform of the quark jet mass distribution 
at hard process scale Q. The result obtained for this function, again to next-to-leading 
logarithmic accuracy, was 
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N f , (2.3) 



and a s defined in the MS scheme. 

In ref. |TJ[], the region of low q 2 in Eq. ( J2.2| ) was ignored, on the grounds that its 
contribution is subleading. We now include it in the following way: we first subtract the 
perturbative contribution from the region q 2 < fi 2 , and then add it back again using the 
non-perturbative effective coupling a e s(q). Changing the order of integration, the amount 
added back is 



6mjim=^r d ^) r-( 

n Jo q Jq 2 /Q 2 u v 



g/Q du 



-uvQ 2 



(2.4) 



Since vQ 2 is conjugate to 1 — T, for 1 — T ^> fii/Q we can safely expand the exponential 
to first order, to obtain 



6 In Jl{Q 2 
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The term involving B(a s ) gives a correction of order 1/Q 2 , which we neglect, together with 



terms of order a 2 s . 



For the perturbative subtraction, we use the next-to- leading-order expansion of a s (q) 
in terms of a s (/x R ), /i R being the chosen renormalization scale: 



a s (q) = a s (/i R ) + ^ In — a 2 s (^ R ) 
Ztt q 



(2.6) 
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where /3 = (HC^ — 2Nj)/3. The expression to be subtracted from the right-hand side of 
Eq. (O) is thus 



which gives 
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Substituting in Eq. (|2.1|) , we see that the leading non-perturbative effect on the thrust 
distribution is simply to shift the perturbative prediction to lower thrust, by an amount 
proportional to l/Q: 

F(T) = F pcrt (T-5T) (2.9) 

where 



5T 



4C F /i : 
Q 
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a (ft) - a s (/x R ) - 



f (ln^ + f + l)«s 2 W • (2-10) 
for the non-perturbative shift in the mean 



This is precisely the formula derived in ref. 
thrust, S (T). 

We remind the reader that the simple prediction (|2.9| ) applies only in the region 1 — T ^> 
Hi/Q, where ft marks the scale below which a e s starts to deviate from a s . For 1 — T ~ ft/Q 
one could explore the effects of retaining more terms in the expansion of the exponential 
function in Eq. Q2.4|) , but this would require a detailed parametrization of a e s, and we have 
not tried it. In practice, we excluded the region 1 — T < 0.05 from all our comparisons 
with data. 



The prediction ( |2.9| ) also does not strictly apply at low values of the thrust, where 
terms involving powers of ln(l — T) are not dominant. In ref. |14]], however, it was found 
that exponentiation has significant effects throughout the region in which the cross sec- 
tion is substantial. In earlier comparisons with data, good agreement was obtained (after 
hadronization corrections) by using a 'logi?' matching scheme, in which essentially all 
known higher-order corrections are exponentiated. It therefore appears natural to adopt 
this scheme and to extend our comparisons with Eq. ( |2.9j ) rather far into the low-thrust re- 
gion. In fact we find good agreement down to T = 0.65, which is even outside the three-jet 
region (T > 2/3). 



For _F pert we take the log i?-matched resummed expression in ref. [14] with no modi- 
fication of the logarithmic terms, i.e. with L = — ln(l — T). Initially, we set the renor- 
malization scale /i R = Q and the infrared matching scale /i R = 2 GeV, as in ref. 0. The 
resulting predictions were compared with data on the thrust distribution in the interval 
0.05 < 1 - T < 0.35 at energies 14 < Q < 161 GeV, as listed in Table 1. 

The best fit values of the two free parameters are 



Ag = 0.235 ± 0.017 GeV , a (2 GeV) = 0.46 ± 0.02 



(2.11) 



(95% confidence level). The corresponding curves are shown in Fig. |l|. The fitted values 
of the parameters are somewhat correlated, as shown in Fig. |2|. The value of A™ corre- 
sponding to a s (Mz) = 0.1186 ± 0.0013, is in good agreement with that obtained by other 



3 



methods. The value of «o is somewhat smaller than, but within two standard deviations 
of, that obtained in ref. §, a (2 GeV) = 0.52 ± 0.03. 

The quality of the overall fit is remarkable (x 2 /d.o.f. = 134/114 = 1.18) - much better 
than those typically obtained^ when hadronization effects are estimated from Monte Carlo 
models ||T6| , [T7| . Furthermore, we see from Table 1 that a large contribution to the x 2 comes 
from the data at 14 GeV, where our fitting region is perhaps too large (i^i/Q = 0.14 > 0.05) 
and there may be complications due to quark mass effects, heavy quark decays and higher 
power corrections. However, the 14 GeV data are valuable because they provide the longest 
lever arm for distinguishing between inverse power and logarithmic energy dependence. If 
we exclude the 14 GeV data altogether, the best fit parameter values do not change, but 
the errors are doubled (dashed curve in Fig. ^). 

To study the dependence on the infrared matching scale, we also performed a fit at 
/ii = 3 GeV. The best fit value of and the quality of the fit did not change significantly, 
and we obtained cio(3 GeV) = 0.374 ± 0.010, again within two standard deviations of the 
value obtained in ref. §, a (3 GeV) = 0.42 ± 0.03. 

We also investigated the dependence on the the renormalization scale /i R , in the range 
Q 2 /2 < n 2 n < 2Q 2 . The best fit parameter values varied from A§ = 0.204 GeV, a (2 GeV) = 

0.457 to = 0.270 GeV, a {2 GeV) = 0.466, respectively. Thus the error in A^ is still 
dominated by the systematic error due to renormalization scale dependence, and our overall 
estimate of this parameter is 

Ag = 0.235 ± 0.035 GeV , a s (M z ) = 0.1185 ± 0.0025 . (2.12) 



Table 1: Data sets and fit results for A^ = 0.235 GeV, a = 0.46. 



Collab. 


Q/GeV Ref. Pts. 


x 2 


TASSO 


14 


18 


8 


28.5 


TASSO 


22 


18| 8 


11.1 


TASSO 


35 


18 


8 


2.1 


TASSO 


44 


18 


8 


7.6 


AMY 


54 


19 


6 


12.5 


OPAL 


91.2 


[20 


30 


19.0 


ALEPH 


91.2 


\>l 


11 


12.8 


DELPHI 


91.2 


£2| 


13 


23.6 


SLD 


91.2 


m 


6 


3.5 


OPAL 


133 


m 


6 


2.9 


DELPHI 


133 


m 


6 


8.4 


OPAL 


161 


m 


6 


2.3 


Total 






116 


134.1 



^Compare, for example, the best fit x 2 /d.o.f. = 2.3 obtained in ref. [^(J at a single energy for the more 
restricted interval 0.06 < 1 — T < 0.30. 
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Figure 1: Data on the thrust distribution at various energies, and the prediction ( |2.9| ) for 
A<g = 0.235 GeV, a = 0.46. The twelve data sets are as listed in Table 1 (from top to 
bottom, multiplied by the factors indicated). 

3 Relation to dispersive approach 

Here we study the non-perturbative contribution to the thrust distribution from the view- 
point of the dispersive approach proposed in ref. ||. As discussed there and in ref. ||, the 
effects of soft gluons on event shape variables are equivalent to a running effective coupling 
only after certain kinematic approximations, which we can clarify using this approach. 

For the kinematics, we take p + and p_ along the initial quark and antiquark directions 
and set 2(p + p_) = 1 (i.e. we measure all momenta in units of Q in this section). The 
light-cone (Sudakov) decomposition of the momentum of an object with mass m is then 

k± + m 2 

k = zp + + ap_ + k± , a = . (3.1) 

Z 

3.1 Exponentiation of "massive" soft gluons 

Each final (massless) soft parton / contributes min{a/, Zf} to (1 — T). This statement is 
based on neglecting quark recoil: kj_r as compared with a/ = k\Jzf. 
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41 < GeV ) 

Figure 2: Best fit values and 95% confidence region for the two fitted parameters. 
Solid/dashed ellipses: including/excluding 14 GeV data. 



The first step consists of assembling final partons into subjets generated by primary 
gluon radiation off the quark-antiquark line and substituting a virtual ( "massive" ) gluon % 
for each subjet. It is implied in doing so that all the secondary partons belonging to a given 
gluon subjet have the same sign of Zf — etf, that is, they lie in the same hemisphere. We 
refer to this as the moving-along assumption. If it is true (say, a/ < z/, i.e. right-movers), 
the total sum = ctj can be attributed to the primary gluon i. (Notice that a 

massive parton contributes min{aj, z^\ to 1 — T as well as a massless one.) This makes 
internal jet structure insignificant and the problem essentially Abelian. 

Interchanging the definitions of z and a for the left-moving virtual gluons, so as to have 
min{«j, Zi} = ctj, we apply the identity 

^i-r-W^^^II^- (3-2) 

i i 

This, together with the factorized matrix element for multiple soft gluon radiation, results 
in a factorized 2-dependence and makes exponentiation straighforward: 

1 da f dv * , m\ t— * / \ i . . 

The radiator function R{y) is obtained from the single soft gluon emission probability 

dz f da Cp 



r dz r da U F r 

J Z J a 7T J 





a s (-m 2 ) 
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(3.4) 



6 



with m 2 > the gluon invariant mass- squared. Here for the sake of simplicity we have kept 
only the double-logarithmic part, corresponding to the main term A(a s ) in ( |2.2| ). 

The crucial moving-along assumption may be correct for sufficiently small gluon angles 
(say, Zi > n 2 a>i with k ~ e) but obviously fails for Zi ~ ai (k\\i — 0, "transversal" gluons). 
We have to assume these do not modify the nature of the leading power correction though 
they may contribute to its magnitude. 

Hereafter we introduce the arbitrary kinematical cut k, to quantify their effect. It 
should disappear after a full treatment of a large-angle gluon emission with decay products 
falling into opposite hemispheres. 



3.2 Running coupling in the radiator 

The exponent formally contains two pieces: the "real" gluon contribution and that from a 
cut virtual one: 



Riy) 



C 



^ /' 1 d m 2 L s (0)5(m 2 ) + ^ 
7T Jo m 2 



1 r 1 dz c z l^ da 


e~ m -1 


J J ran Z Jm 2 /z Ct 





(3.5) 



The lower limit a > m 2 /z comes from k\ > 0. The upper limit a < z/k 2 is the moving- 
along condition enhanced by the factor k > 1 as discussed above, and so is the lower limit 
in z: z/k 2 > a > m 2 /z. It is the latter that brings in non-analyticity in m 2 , which is 
crucial for generating power corrections [@,[|. 

To simplify the analysis we may differentiate with respect to v and perform the a 
integration: 



vRJv) 



C 



1 /I dm 2 {a s (0)5(m 2 ) + P(m2) 
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1 dz r 

ran Z 



The next step is to integrate with respect to m 2 by parts using (see ref. |J) 

dm 2 

This gives 



2 pO 2 ) 



da e ff(m 2 ) . 



(a,(0) - aeff(O)) 
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1 ( v\ dz 



-vm 2 1 z 



Z J Z 



The first term vanishes identically since 

dm 2 



a s (k 



o m? + k 2 

poo dfn 2 f°° 
a s (0) = - — 7rp{m 2 ) = - da eS (m 2 ) = ct s (oo) - a eff (0) = a cfT (0) 
Jo m z Jo 

Performing the ^-integration we finally arrive at 



vRJu) 



C F f^dm 2 
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-otesvm 
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(3.7) 



(3.8) 



(3.9) 



When evaluating (|3.9|) perturbatively, we do not distinguish between a eff and a s = a e ff(l + 
O(otg)). This results in the above expression ( j2.2| ) with the leading term A(a s ) in the 
physical scheme of refs. [|Tl~| , |2"7| for the coupling. The next-to- leading term B(a s ) is also 
easy to reproduce by keeping non-soft contributions in the elementary radiation probability 
(0). 

In the non-perturbative region m 2 < p 2 we can expand the exponentials to obtain the 
result in Eq. (|2.9|) with the shift expressed in terms of the first non-analytic moment of 
the effective coupling (A± in the notation of ref. ||). The only difference from Eq. (|2.5| ) is 
an overall factor of 1/k, indicating that the magnitude of the shift is sensitive to the soft, 
large-angle region of gluon emission. 



4 Discussion 



We see from the above discussion that the non-analyticity in m 2 in Eq. ( |3.9|) can be traced 
back to the very kinematical region z ~ a ~ m which does not respect the crucial moving- 
along assumption. This does not affect the nature of the leading power correction but 
it could affect its magnitude, since reducing the contribution of this region scales down 
the shift 5T by the factor k which we have introduced to quantify sensitivity to large- 
angle gluon radiation. In the large-angle region one can expect an essential modification 
of the inclusive spectral density p due to the specific kinematics of the thrust: the "decay" 
products of a timelike virtual gluon in this region may make different contributions to the 
event shape, depending on the kinematics of the decay. In the large- Nf model studies of 
ref. ||, this effect was not found to be large for the thrust, but it does depend on the shape 
variable involved. 

In the case of other jet shape variables which have the property of exponentiation, such 
as the heavy jet mass [28| and jet broadening p9| , we expect 1/Q corrections to be generated 



by the same mechanism. However, we found that for these quantities the leading non- 
perturbative effect is not well represented by a simple shift in the distribution. Furthermore, 
the modifications due to the large-angle region discussed above will be different for different 
jet shape observables. Thus for an extension to other related jet observables a quantitative 
analysis of the large-angle region has to be pursued. 



Acknowledgements 

We have benefited from many valuable conversations on this topic with S. Catani, M. 
Dasgupta, G. Marchesini, A.H. Mueller, P. Nason, G. Salam and M.H. Seymour. Yu.L.D. 
thanks the Cavendish Laboratory and B.R.W. thanks the CERN Theory Division and the 
St Petersburg Nuclear Physics Institute for hospitality while part of this work was carried 
out. 



8 



References 



1. For reviews and classic references see: 

V.I. Zakharov, Nucl. Phys. B385 (1992) 452; 

A.H. Mueller, in QCD 20 Years Later, vol. 1 (World Scientific, Singapore, 1993). 

2. M. Beneke, V.M. Braun and VI. Zakharov, Phys. Rev. Lett. 73 (1994) 3058; 
P. Ball, M. Beneke and V.M. Braun, Nucl. Phys. B452 (1995) 563; 

M. Beneke and V.M. Braun, Nucl. Phys. B454 (1995) 253. 



3. G. Grunberg, Phys. Lett. 372B (1996) 121, CPTH-PC463-0896 ||hep-ph/9608T75[ 
N.J. Watson, CPT-96-P-3347 ||hep-ph/ 9606381] 



4. B.R. Webber, Phys. Lett. 339B (1994) 148; see also Proc. Summer School on Hadronic 
Aspects of Collider Physics, Zuoz, Switzerland, 1994 ||hep-ph/ 941 1384 . 



5. Yu.L. Dokshitzer and B.R. Webber, Phys. Lett. 352B (1995) 451. 

6. Yu.L. Dokshitzer, G. Marchesini and B.R. Webber, Nucl. Phys. B469 (1996) 93. 

7. G.P Korchemsky and G. Sterman, Nucl. Phys. B437 (1995) 415; see also Proc. 30th 
Rencontres de Moriond, Meribel-les-Allues, France, 1995 ||hep-ph / 950539 1 . 



8. P. Nason and M.H. Seymour, Nucl. Phys. B454 (1995) 291. 

9. R. Akhoury and V.I. Zakharov, Phys. Lett. 357B (1995) 646, Nucl. Phys. B465 (1996) 
295. 



10. V.M. Braun, NORDITA-96-65P ||hep-ph/96102H 



M. Beneke, SLAC-PUB-7277 [ |iep-ph/9609"2l"5| 



11. Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan, Phys. Rev. D53 (1996) 89. 

12. DELPHI Collaboration, P. Abreu et al, Z. Phys. C73 (1997) 229. 

13. S. Catani, L. Trentadue, G. Turnock and B.R. Webber, Phys. Lett. 263B (1991) 491. 

14. S. Catani, L. Trentadue, G. Turnock and B.R. Webber, Nucl. Phys. B407 (1993) 3. 

15. E. Farhi, Phys. Rev. Lett. 39 (1977) 1587. 

16. T. Sjostrand, Comp. Phys. Commun. 39 (1984) 347; 

M. Bengtsson and T. Sjostrand, Comp. Phys. Commun. 43 (1987) 367. 

17. G. Marchesini, B.R. Webber, G. Abbiendi, I.G. Knowles, M.H. Seymour and L. 
Stanco, Comp. Phys. Commun. 67 (1992) 465. 

18. TASSO Collaboration, W. Braunschweig et al., Z. Phys. C47 (1990) 187. 

19. AMY Collaboration, Y.K. Li et al., Phys. Rev. D41 (1990) 2675. 

20. OPAL Collaboration, PD. Acton et al., Z. Phys. C59 (1993) 1. 

21. ALEPH Collaboration, R. Barate et al, CERN-PPE-96-186 (1996). 

9 



22. DELPHI Collaboration, P. Abreu et al., CERN-PPE-96-120 (1996). 

23. SLD Collaboration, K. Abe et al, Phys. Rev. D51 (1995) 962. 

24. OPAL Collaboration, G. Alexander et al, Z. Phys. C72 (1996) 191. 

25. DELPHI Collaboration, P. Abreu et al., CERN-PPE-96-130 (1996). 

26. OPAL Collaboration, K. Akerstaff et al., CERN-PPE-97-015 (1997). 

27. S. Catani, G. Marchesini and B.R. Webber, Nucl. Phys. B349 (1991) 635. 

28. S. Catani, G. Turnock and B.R. Webber, Phys. Lett. 272B (1991) 368. 

29. S. Catani, G. Turnock and B.R. Webber, Phys. Lett. 295B (1992) 269. 



10 



